function shares = housing_supply(rs, Hs)
    global elasts
   
    shares =  ( Hs(1:end) .*  rs(:)  ).^ elasts(1:end);
 
end